Normal state of a polarized Fermi gas at unitarity 
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We study the Fermi gas at unitarity and at T = by assuming that, at high polarizations, it is 
a normal Fermi hquid composed of weakly interacting quasiparticles associated with the minority 
spin atoms. With a quantum Monte Carlo approach we calculate their effective mass and binding 
energy, as well as the full equation of state of the normal phase as a function of the concentration 
X — ni/ni of minority atoms. We predict a first order phase transition from normal to superfluid 
at Xc — 0.44 corresponding, in the presence of harmonic trapping, to a critical polarization Pc = 
{Nf — Ni)/{Ni + Ni) — 77%. We calculate the radii and the density profiles in the trap and predict 
that the frequency of the spin dipole mode will be increased by a factor of 1.23 due to interactions. 
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Recent experiments on degenerate gases of ^Li with 
a mixture of two hyperfine species have explored the 



we shall assume that the change in E is given by 



and have led to a num- 
1 . One of the ma- 



physics of Fermi gases 0, 
ber of theoretical analyses 
jor experimental observations has been the occurrence of 
phase separation if the mixture contains more atoms of 
one species than of the other, i.e., if the gas is polarized. 
Some of the experiments suggest that in the unitary limit 
of strong interactions there are three phases: an unpo- 
larized superfluid phase, a mixed phase which exhibits a 
partial polarization and a fully polarized gas. We now 
have a good understanding of the superfluid phase which 
has been the subject of numerous theoretical and experi- 
mental studies while the fully polarized phase is an ideal 
Fermi gas since atoms in the same spin state do not in- 
teract with each other. However, for intermediate polar- 
izations, when both species are present, the nature of the 
mixed phase is not understood. 

Here we study the mixed phase in the unitary limit 
by adopting an approach inspired by the theory of dilute 
solutions of '^He in ''He '3]. We will assume that the 
majority species (j) forms a background experienced by 
the minority species (j) and that the latter behaves as 
a gas of weakly interacting fermionic quasiparticles even 
though the t ^ i atomic interaction is very strong, being 
characterized by an inflnite scattering length. In other 
words, we will assume that the system is a normal Fermi 
liquid, which will allow us to characterize the energy of 
the gas in terms of a few parameters and, by calculating 
these with a quantum Monte Carlo approach, allow us to 
make various predictions of experimental relevance. 

We begin by writing the expression for the energy E of 
a homogeneous system in the limit of very dilute mixtures 
and at zero temperature. The concentration of J, atoms 
is given by the ratio of the densities x = ni/n^ and we 
will take it to be small. If only f atoms are present then 
the energy is that of an ideal Fermi gas E{x = 0) = 
3/5Ef^N^ , where iV| is the total number of t atoms and 
Epi — /i^/2m(67r^n|)^/'^ is the ideal gas Fermi energy. 
When we add a | atom with a momentum p (|p| <C pft)j 
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In other words, the | atom in the t gas behaves as a 
quasiparticle with a quadratic dispersion and an effec- 
tive mass m* . In addition, there is a "binding" energy 
—3/5Ep^A of the | atom to the Fermi gas of t atoms. 
This binding energy must be proportional to Ep^ since 
there is no other energy scale in the unitary limit and we 
have used the factor 3/5 for later convenience. We shall 
further assume that this quasiparticle is a fcrmion 

When we add more | atoms, creating a small finite den- 
sity ni, they will form a degenerate gas of quasiparticles 
at zero temperature occupying all the states with mo- 
mentum up to the Fermi momentum pp^ = fi(67r^n|)^/^. 
The energy of the system can then be written in a useful 
form in terms of the concentration x as: 
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Eq.(|2Il is valid for small values of the concentration x, 
i.e. when interactions between | quasiparticles as well as 
further renormalization effects of the parameters can be 
neglected. 

In the following we will calculate A and m* using a 
fixed-node diffusion Monte Carlo (FN-DMC) approach, 
already employed in earlier studies 0, 0| . We use the 
same attractive square-well potential to model the inter- 
actions between [ and f atoms: V{r) — — Vq for r < Rq 
and V{r) = otherwise. The short range Rq is chosen 
as 2n|i?g = 10~^. The depth Vq is fixed by the unitar- 
ity condition |a| = oo for the s-wave scattering length a 
and corresponds to the threshold for the first two-body 
bound state in the well: Vq = tt^/i^/ (4mi?o). For a single 
I atom in a homogeneous Fermi sea of t atoms the trial 
wave function ipT, which determines the nodal surface 
used as an ansatz in the FN-DMC calculation, is chosen 
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to be of the form [15 



using the trial wave function 
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where denotes the position of the | atom and r^^ = 
|rj^ — r^i I . In this equation the plane wave exp (iq • ) cor- 
responds to the impurity travelling through the medium 
with momentum fi,q = {nj.,ny,nz)27Th/ L, where L is the 
lenght of the cubic box and the nt are integers describing 
the momentum in each coordinate. Furthermore, £'|(iV|) 
is the Slater determinant of plane waves describing the 
Fermi sea of the Nf atoms and the Jastrow term /(r) 
accounts for correlations between the impurity and the 
Fermi sea. The correlation function /(r) is chosen as 
in Ref. 01 ■ We consider a system of = 33 atoms 
and periodic boundary conditions. From the energy E 
of the system with N-^ + 1 atoms we subtract the exact 
energy Eo{N^) of the Fermi sea. The result is shown 
in Fig. Has a function of {q/kp])^- From a linear best 
fit we obtain the following values: A = 0.97(2) and 
m* /m — 1.04(3) which is consistent with the general 
inequality m* /m > 1. We have checked that finite-size 
corrections associated with the number N-^ of atoms in 
the Fermi sea are below the reported statistical error. 
It is worth noticing that the binding energy of a single 
I atom almost coincides with the average energy of an 
atom in the Fermi sea. A similar coincidence was found 
for the superfluid gap A at unitarity [l^ . 

A relevant question is to understand whether the equa- 
tion of state (0) is adequate to describe regimes of large 
values of x where interaction between quasiparticles and 
other effects might become important. To answer this 
question we have carried out a FN-DMC calculation of 
the equation of state at finite concentrations x ~ Ni/N-^ 




FIG. 1: Excitation spectrum of a | atom in a Fermi sea of | 
atoms. The dashed line is a linear best fit to the FN-DMC 
results from which we extract the values A and m*. 
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where i and i' label, respectively, t and | atoms. The 
nodal surface of the wave function tpT is determined by 
the product of Slater determinants D-^{N-i)Di{Ni) and 
coincides with the nodal surface of a two-component ideal 
Fermi gas. As a consequence, the wave function in Eq. 
Q is incompatible with off-diagonal long-range order 
(ODLRO) and describes a normal Fermi gas. In contrast, 
the BCS-type wave function used in Refs. 0,0] is com- 
patible with ODLRO and describes a superfluid state. 
A direct comparison between the groundstate energy of 
the normal and superfluid states can be carried out for 
equal numbers of t and | atoms, = N^, with the re- 
sult Esf/{3/5Ef^N^) 0.84(2) and £;n/(3/5£;ft^^t) = 
1.12(2) showing the instability of the normal state for 
X — 1 (see Fig.|2J. 

The results for the equation of state of the normal 
Fermi gas are shown in Fig. [2 To reduce finite-size 
effects we have considered closed-shell configurations 
Arx=7,19,27,33 with iV|=27,33. In Fig. |21we also show 
the prediction of Eq.(|2l based on noninteracting quasi- 
particles (dashed line). For small values of x we find very 
good agreement, but for larger concentrations effects of 
interactions between quasiparticles start to be important 
and deviations from Eq.Q become visible. The solid line 
is obtained from a polynomial best fit to the FN-DMC 
results. 

From the equation of state of the mixed phase we can 
determine the transition between the fully polarized and 
the mixed phases as well as the transition between the 
mixed and the unpolarized superfluid phases The 
equilibrium condition is obtained by imposing that the 
chemical potential and the pressure be the same in the 
two phases. It is useful to express the results in terms of 
the chemical potential = (^|+/i4)/2 and of the effective 
magnetic field h = — ni)/2. Here — dE/dN^i 
is the chemical potential of each spin species. The tran- 
sition between the fully polarized and the mixed phase 
is second order and takes place at a; = where we find 
Mi/Mt = —3/5^, corresponding to /i//Lt = 3.78. The tran- 
sition between the mixed and the unpolarized superfluid 
phases is instead first order and is simply obtained via 
the standard Maxwell construction considering the tan- 
gent to the equation of state of Fig. |21 crossing the super- 
fluid point EsyK'^/^Ef^N^) = 0.84 at x=\. We obtain 
the critical value Xc = 0.44 at the transition. For smaller 
values of x the system remains in the normal state, while 
above the critical concentration Xc the system will begin 
nucleating the superfluid and phase separate into those 
two states. The phase transition is characterized by the 
value /ix//^t = 0.017 corresponding to h/^ = 0.96 [l5| . 
Note that at the critical value Xr = 0.44 the difference 
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between the best fit and the prediction of Eq. ^ is quite 
small so that this latter energy functional describes well 
the whole normal phase. This supports the idea that the 
normal state is a gas of weakly interacting quasiparticles. 
Eq.® would actually predict the value Xc ~ 0.50 for the 
critical concentration, which corresponds to h/fi = 0.93. 

It is useful to compare the above results with BCS 
theory which also predicts 3 phases at unitarity 0, 
a superfluid with energy _EBCs/(3/5i?_F|iV|) = 1.18, a 
mixed state which is a noninteracting partially polarized 
gas and a fully polarized gas. The mixed state energy is 
simply the kinetic energy and is an increasing function 
of X. The tangent to the curve is at Xc = 0.04, cor- 
responding to fJ']'/fJ'i ~ 0.1, thereby leading to a much 
reduced normal region with respect to the predictions of 
our FN-DMC calculation. 

Now we turn to the trapped case. While we will dis- 
cuss the situation only at zero temperature we should 
note that temperature can have an important effect on 
the density profile of the | atoms. In experiments we 
usually have the condition ksT ^ Ep^. But since 
n|, we might be in a situation where ksT > Epi 
and would therefore need to use the appropriate thermal 
distribution. With this caveat we turn to the inhomo- 
geneous situation where we shall use the local density 
approximation (LDA) [l^ . In a harmonic trap with po- 
tential V{r) = X^i "^^i the local chemical poten- 
tials become /^-fj.(r) = fi^^ — T^(r). For small concen- 
trations (in particular x <^ Xc), where only the normal 
state is present and where we can neglect the change 
in /x-f due to the attraction of the | atoms, n-f is the 
Thomas-Fermi density of an ideal gas whereas is a 




FIG. 2: Equation of state of a normal Fermi gas as a function 
of the concentration x (circles) . The solid line is a polynomial 
best fit to the FN-DMC results. The dashed line corresponds 
to Eq.l|2l. The dot-dashed line is the coexistence line between 
the normal and the unpolarized superfluid states and the ar- 
row indicates the critical concentration Xc above which the 
system phase separates. For x = 1, both the energy of the 
normal and of the superfluid (diamond) states are shown. 
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FIG. 3; Radii of the three phases in the trap in units of the 
radius i?| = aho(48A''|)^''^ of a noninteracting fully polarized 
gas, where aho is the harmonic oscillator length. 



Thomas-Fermi profile with a modified harmonic poten- 
tial V ^ V X {1 + 3/5 A). The potential seen by the 
I atoms is more confining due to the attraction to the 
t atoms. This also has consequences for the collective 
modes of the system: it leads us to predict that the spin 
dipole mode - the mode where the | atoms oscillate as 
a whole in the midst of the | atom cloud - will have a 
frequency given by 
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So, a direct measurement of the oscillation frequency of 
the I atoms in a dilute mixture with a; <C 1 would provide 
a useful test of this Fermi liquid theory and in particular 
of the numerical estimate of the parameters A and m* 
p7j . For larger concentrations the system will exhibit 
a central superfiuid core with chemical potential hsf = 
(/Z| + ijLi)/2, whose radius in a spherical trap is given by 
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where is the radius of the j component and ^J^i/ ^jl] is 
calculated at the transition point. In Fig. El we plot the 
radius Rsf of the superfluid component together with 
the radii of the minority and majority components, i?| 
and i?| respectively, in units of R^ — aijo(487V|)^/^ as a 
fimction of the polarization P = {N^ ~ N{)/{N^ + N{) 
of the sample. We predict that the superfluid phase 
disappears at Pc = 0.77 in good agreement with the 
experimental findings of {3|. Notice that as P 1, 
i?X ~ i?0[(l + 3/5yl)m7m]-'/^[(l-P)/(l + P)]'/^ ^ 
while i?| approaches the noninteracting value. In the 
opposite P — !■ limit the radii converge to the known 
value {fi/ Epi^y/'^ ~ 0.80. It is worth noticing that the 
BCS approach would predict the value Pc > 0.99 at uni- 
tarity 5, 6] pointing out the dramatic role played by the 
binding energy in the mixed state. 



4 



0.8 




r 



FIG. 4: Density profiles (soUd lines) for the two spin com- 
ponents in a spherical harmonic trap at the critical total po- 
larization Pc — 0.77. The density profiles (dashed lines) for 
a noninteracting Fermi gas with the same polarization are 
also shown. The inset shows the density difference profile. 
Densities are given in units of the central density of the non- 
interacting gas. The radial coordinate is given in units of iij 
(see Fig. . 



At the transition between the superfluid and the mixed 
normal phase the density of the two spin species ex- 
hibits a discontinuity, revealing the first order nature of 
the transition. The densities n| and jump from the 
superfluid value hsf to the values ~ l.Olnsp and 
ni = XcTi]' ~ OAAnsF respectively, as one enters the nor- 
mal phase. The discontinuity is an artifact of LDA and 
the inclusion of surface energy effects could In Fig. 0] we 
plot n-\ and ni as well as the difference n-^ —ni as a func- 
tion of position in a spherical trap when P — Pc. These 
results, based on LDA, apply also to anisotropic traps 
through a simple scaling transformation. We find that 
both the total density n-f -I- and the density difference 
ri| — ni increase monotonically towards the center [l8l| . 
If P < Pc then n-f = in the central superfluid region. 

In conclusion, we study the polarized Fermi gas at uni- 
tarity as a normal Fermi liquid composed of weakly inter- 
acting quasiparticles associated with the minority atoms. 
These have a quadratic dispersion and, in the x — * 
limit, have an effective mass m* ~ 1.04(3)m with a bind- 
ing energy -3/5£'ft x 0.97(2), calculated with a FN- 
DMC approach. We derive an energy functional Eq.Q 
with those parameters, assuming noninteracting quasi- 
particles and, using FN-DMC calculations at higher val- 
ues of X, show that corrections to the energy functional 
are small, even for relatively high concentrations. As- 
suming that no polarized superfluid phases exist, we pre- 
dict a normal/superfluid first order phase transition at 
a critical value Xc — 0.44, corresponding, in the pres- 
ence of harmonic trapping, to a critical total polarization 
Pc = 0.77. We also predict for small concentrations an 
increase by a factor of 1.23 of the frequency of the spin 
dipole mode with respect to the noninteracting value. A 



further application of the equation of state of the mixed 
phase could be, for example, the calculation of the dipole 
polarizabihty of the trapped Fermi gas [l^ . 
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